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ABSTRACT 

It is shown that major discrepancies exist between experiments and theory 
for ice crystal growth from solution. Accurate data, taken in a microgravity 
environment, approximate analytical models, and exact (probably numerical) 
models all are needed to advance our understanding of ice crystal growth 
phenomena. 

A new approximate semi -empirical theory is presented which predicts that 
a relatively sharp transition from natural convection control to diffusion 
control for ice growth in pure water occurs at a subcooling of about 10°C (a 
reduced temperature difference of about 0.125). No reliable data exist to 
test this prediction. The theory also predicts qualitatively the growth of 
ice in NaCl solution in which maxima in the growth rates are observed at 
various levels of subcooling. 

For pure water the exponent on the subcooling (n) predicted by 
Equation (11) when cast in the form of Equation (1) is 1.9 , which is in 
reasonable agreement with the experimental value. 
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The purpose of the proposed single crystal ice growth experiment is to 
clarify the roles of thermal and concentration driven convection in crystal 
growth by measuring growth rates, crystal shape, and temperature and 
concentration fields in the neighborhood of the single crystal, in the 
microgravity environment of space, by using shadowgraphs, Schlieren and 
other optical techniques. It is expected that these data, taken in certain 
critical ranges of the operating parameters, will help to explain the more 
serious discrepancies which exist among the various theories of crystal 
growth and the best available data. 

Ice crystal growth rates usually are represented by an equation of 
the form 

V = aAT" (1) 

where v is the growth velocity and aT is the subcooling. Ivantsov's theory 
predicts n = 3 for parabolic platelets. The best experimental values of 
n for growth in quiescent water are closer to 2.17. Recent experimental 
results of Barduhn and Kallungal (1977) have shown ice growth rates to 
depend on Crystal orientation with respect to the gravitational field, and 
observed rates are an order of magnitude larger than Ivantsov's theory 
predicts. 

Growth rates in the a-axis direction are controlled by heat and mass 
transfer; crystal formation kinetics play no part. Observed growth rates 
depend on the salt concentration of the solution and growth rate maxima are 
observed with about 1% salt solutions. All of this provides convincing 
evidence that thermal and concentration driven convection play a major role 
in the growth of ice crystals despite their small size. However, no adequate 
theory exists to describe these phenomena. 



; There are fundamental differences between the best data available on 
ice and succinonitrile and the theories of crystal growth which have been 
developed. The data of Barduhn and Kallungal (1977) for ice and those of 
Glicksman, Shaefer and Ayers (1976) for succinonitrile show that the tip 
radius, R, varies as the reciprocal of under cooling so that 

R = AaT^ (2) 

where m = -1. However the linear stability theory of Langer and Muller- 
Krumbhaar (1978) predicts a much larger negative value for m which is 
closer to -2. This is shown in Langer and Muller-Krumbhaar's Figure 4 which, 
being a log-log plot, tends to make the differences between theory and 
experiment less obvious. On the other hand the maximum velocity principle 
predicts m = -1 , which is in accord with experimental data. Unfortunately 
the maximum velocity assumption significantly underestimates A. 

The Langer- Krumbhaar theory predicts the tip velocity data of Glickman 
et al . better than the tip radius data with deviations tending to get larger 
as under cooling increases. This seems to illustrate that linear stability 

theory often is a useful guide, but seldom is a definitive statement on the 

> 

quantitative behavior of real systems. 

Fernandez and Barduhn (1967), Vlahakis and Barduhn (1974), Simpson 
et al. (1974), Kallungal and Barduhn (1977) and Huang (1975) have studied 
various aspects of ice crystal growth in forced convection systems with 
velocities up to 70 cm/sec. Equation (1) continues to apply here except 
n = 3/2 as predicted by Fernandez using boundary layer theory (and the 
maximum velocity principle) and observed experimentally provided the forced 
velocity is large enough. Equation (2) also applies to forced convection; 
as before m is observed and predicted to be -1. However the maximum 
velocity principle again predicts values of A which are too low by a 
factor of 3 or 4 for ice. The data show conclusively that tip radius is 
independent of forced flow velocity over six orders of magnitude; it also 
is independent of crystal orientation. 

The data on ice growth in a quiescent ambient fluid deafly are 
influenced by natural convection as evidenced by their sensitivity to 
crystal orientation. No adequate theory, which includes thermal . convection, 
has ever been proposed to describe ice crystal growth even in pure water 
and the situation is complicated markedly when the crystals are grown in 
saline solution. The difficulty of constructing an exact theory, in the 
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sense of obtaining exact solutions of the transport equations, is illustrated 
by the interesting article of Kroeger and Ostrach (1974) on a model which 
simulates a casting problem and includes forced and natural convection. 

Models which include exact solutions of the relevant transport equations 
are needed and should be worked on. However, in the meantime it may be 
possible to gain further insight into crystal growth, including thermal 
convection, by constructing approximate models which perhaps are viewed best 
as "order of magnitude estimates" of the real phenomena occurring. One such 
possible approximate approach, which appears to be new, and which includes 
both thermal convection and- the moving boundary effect will be sketched here. 

Consider for illustration a crystal growing as a parabolic platelet. 

A heat balance on the crystal, which will be assumed to be. isothermal , gives 




Nu 


k(T. - T ) 

' T 00 ' 


RP^ L 



UR 

where Nu = Nusselt number (-jp). k = thermal conductivity of the solution, 

L = latent heat of fusion, = crystal density, T.. and = interfacial and 
ambient temperatures of the fluid. From the Ivantsov theory, as discussed 
by Horvay and Cahn (1961), if p < < 1, we find 

where Cp = heat capacity of the melt and a = thermal diffusivity. Thus, if 
we assume density differences between the crystal and solution are small. 
Equations (3) and (4) lead to 

" %.B. " I L/Cp 

where the subscript M.B. denotes moving boundary and this Nusselt number 
reflects two dimensional heat conduction in the melt with the "apparent 
convection" created by a moving boundary. 

Obviously (5) tells only part of the story; it neglects thermal 
convection entirely. For a parabolic platelet, even with a fixed boundary, 
it is a formidable job to solve the thermal convection problem with the 
two dimensional conduction which must be included in the low Rayleigh 
number (Ra) systems encountered in crystal growth. However, some information 
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can be obtained very easily from the experimental results of Gebhart and 
Pera (1970) on heat transfer at low Ra from circular cylinders of very small 
diameter. These data should be relevant to ice growth if heat transfer at 
the stagnation point of the tip of the parabolic platelet does indeed control 
the growth process. The correspondence between Fernandez's forced convection 
theory, which employs the tip radius as the characteristic dimension, and the 
experimental growth rate data lends considerable credence to this assumption. 
One can show (for Ra < 10'^) from Gebhart and Pera's data that the Nusselt 
number due to two dimensional conduction and natural convection is 


Nu 


N'.C. 


0.47 Ra 


0.075 


( 6 ) 


where 


Ra = 


egR" (Ti “ TJ 


\>a 


e = thermal conductivity, g = gravitational constant, v = kinematic viscosity. 

Equation (6) is very different from the conventional boundary layer result, 
1/4 

Nu^j c ~ R® > which grossly underestimates heat transfer rates at very 
low Ra. If one uses the boundary layer approach, then an arbitrarily enlarged 
length scale is needed to obtain even qualitative agreement between theory and 
experiment. 

Churchill (1977) in numerous papers has discussed the accuracy of combining 
rules for finding the Nusselt numbers for "mixed mode convection!' among other 
things. Ruckenstein (1978) has provided a theoretical justification for 
Churchill's procedure. If we consider natural convection and the moving 
boundary as two processes competing for dominance, which seems plausible but 
doesn't seem to have been done before, Churchill's work would suggest that a 
reasonable approximation can be obtained from 

Nu3 = 


or 


Nu = [n4b f 

Therefore, one obtains using (3),' (5), (6), and (7), 




4 L/C 
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(7) 


( 8 ) 


where at = T^ 


Obviously (8) cannot be viewed at this stage as anything 
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more than a plausible speculation. However, Equation (8), and other 
expressions like it which can be derived in a similar fashion, deserve the 
effort of being compared with experimental data and perhaps modified 
accordingly. 

It is worth noting that if NUu « = 2Nu„ - , then Nu„ « “ 0.96Nu. 

n.D. ri.u- n.D* 

Even if the factor is 1.25 rather than 2, one gets NUn^j g = 0.87Nu. That is, 
when either of the Nusselt numbers is larger than the other, the mode 
representing the larger one dominates the process. This can be used to 
determine the aT above which natural convection is negligible by setting 
Nu^^ g ^ C solving for aT to obtain aT = 10®C. If we had used 
boundary layer theory, we would have obtained aT = 0.5°C which, on the basis 
of Kallungal's data, is too small. Unfortunately, reliable data for aT, 
equal to 10 or more, do not exist to check this prediction. It also is worth 
noting that Nu„ ^ , at the very low Ra, is insensitive to the value of Ra. 

N .0 . 4 

For example, Nuj^ ^ changes only by a factor of 2, while Ra changes by 10 . 
However, depending on ones purposes, this change may be enough, as noted 
above, to cause Nu to be dominated by NUj^ g . 

The forced convection ice crystal growth data for pure water are in good 
agreement with the theory of Fernandez and Barduhn (1967) if the flow 
Reynolds number based on tip radius is high enough. However, the experimental 
values of tip radius are about four times larger than the theory predicts. 

That is, theory and experiment show that m, in Equation (2), equals minus 
unity, but the experimental value of A = 0.6, is about four times the 
theoretical one; the reasons for this are unknown. 

In the subcooling range of 0°C to 1°C, Kallungal and Barduhn (1977) 
found experimentally that the forced convection growth theory predicts growth 
rates in pure water down to a Reynolds number, based on tip radius, of 
5 X 10”^. If one examines carefully the data of Gebhart and Pera (1970) 
on mixed convection heat transfer to fine wire circular cylinders, one finds 
that natural convection alters forced convection by about 3% when 

Re = 0.5 6r^/^ 

The Grashof number, based on experimental values of ice crystal tip radius, is 
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Thus, for 0.1 £ AT £ 1, one gets 

1.3 X 10"^ < Re < 4 X 10"^ 


Consequently, on the basis of mixed convection heat transfer data, one indeed 
would expect natural convection to influence forced convection ice crystal 
growth rates in the Re range where the forced convection theory was found to 
fail- 


The problem of growing ice crystals in saline (NaCl) solution is much 
more complicated and no theories have been proposed which even describe the 
experimental data qualitatively. The important facts are: growth rates 

increase substantially up to salt concentrations of about 1%; the tip radiu*' 
no longer is inversely proportional to aT; the location of the density 
maximum depends on concentration. 

Actually, some simple calculations can shed considera 
influence of solute concentration on crystal growth rates, 
together with its counterpart for mass diffusion, yields 


ile li 
Equ? 
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(!k) /Scx 


(T. - T ) 
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( 9 ) 


k R V 

where Sh = Sherwood number (-^)5 Sc = Schmi ' liniber (-q), D = diffusion 
coefficient, = mass transfer coefficien+ ^oo “ interfacial and 

ambient concentrations of solute. Equa^^ (9) shows that and that 

this "concentration polarization^ inr/ ^ses with subcooling. Obviously then, 
the extent of freezing point ceprr^ ion increases with subcooling. If^forced 
convection dominates, to a reas^/^ble approximation (Nu/Sh) « (Pr/Sc) 
where Pr 13.5, Sc = 2600 f/ ice-water-NaGl systems. Since T^ is the 
equilibrium freezing tempera-ire, corresponding to (taking into account 
the Gibbs-Thompson curvaty^ effect and freezing point depression), Equation 
(9) completely determi^ies once T^ and are specified. Furthermore, 
the magnitude of increases, as does the freezing point depression, as 
increases. Thus, with pure forced convection, the only effect of the 
concentration field is to create a freezing point depression, and one expects 
the growth rate for a given subcooling to decrease as increases (subcool ing, 
at is defined as - T , where Tn is the freezing point at C^). In effect, 
as C increases a greater fraction of aT is consumed by the freezing point 
depression and Gibbs-Thompson curvature and less is left to drive the 
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crystallization process. This is precisely what Huang (1975) has observed 
experimentally, and the forced flow velocity necessary for forced convection 
to dominate was between 10 and 60 cm/sec, which is much higher than those 
necessary in the absence of salt. Below this velocity, maxima were observed 
in the plots of growth rates versus C^. 

Clearly, the maxima in the growth rates must be due to natural convection 
induced by concentration differences. Much larger density differences are 
created by concentration than temperature differences and concentration 
differences increase with for a given subcooling. Also, the Rayleigh 
number for mass diffusion is about 200 times that for heat transfer if they 
both are based on the same density difference. All of this means that natural 
convection created by concentration differences will be very strong and will 
increase with C^. Consequently, a competition is set up between the freezing 
point depression effect, which decreases growth rates, and the augmented 
transfer by natural convection, which increases growth rates. Apparently/ 
natural convection wins at C^ less than 1% and freezing point depression 
prevails about 1% salt by weight. 

In the low aT region, natural convection dominates and this is augmented 
by concentration differences. To describe this mathematically, first note 
that 

T, - T. . (Tj - TJ - [Kf (W - W ) + (10) 

where = 31.8 is the freezing point depression constant, and are 
the mass fractions of solute at the interface and in the bulk^ solution, and 
y is the interfacial tension. 

Equations (9) and (10) enable one to calculate which then can be 
used to determi ne the growth vel dci ty from ^ 

W5 

v=^Sh[l-^l (IT) 

^i 

To use Equations (9), (10), and (11), one must know Nu, Sh, and have a 
relationship between R and AT. Realistic theoretical expressions for Nu and 
Sh are badly needed, but it is a formidable job to obtain them. One can get 
a rough estimate for these expressions by assuming that the analogy between 
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heat and mass transfer holds and that Equation (6) applies with the Grashof 
number given by 

Gr = (- V '’* ) (12) 

Pi - Poo 

where — ^ reflects both temperature and concentration differences. This 

leads to “ (Nu/Sh) = (Pr/Sc)^'^^^. Thus, all necessary parameters are known 
once T^ and are specified. 

If one uses (11) with these crude assumptions and the experimental 
expression for R in terms of aT, the results when compared with experimental 
growth data are encouraging. Equation (11) predicts maxima in the growth 
versus curves, as are observed experimentally. Also, the magnitudes of 
V predicted are reasonably close to the experimental values. However, the 
theoretical maxima are not as pronounced as the experimental ones. Predicted 
rates are somewhat too high near = 0 and somewhat too low at concentrations 
between 1 and 2 weight percent. When one casts (11) in the form of (1), one 
finds n = 1.9 for = 0 and n = 2.1 for W^ = 0.02, which are in 
reasonable agreement wTth the experimental values. 

It is not clear why m =-1.5 in Equation (2) for the tip radius as 
observed by Vlahakis and Barduhn (1974) when growth is in a 0.15% NaCl . 
Hydrodynamic effects do not seem likely to account for this because in pure 
water m = -1 over a change of 10-fold in aT and six orders of magnitude 
for the forced velocity of pure water. This question needs further study. 

It is clear that some very important aspects of crystal growth are 
clouded by thermal and concentration driven convection. Accurate data, 
taken in a microgravity environment, approximate analytical models, and 
exact (probably numerical) models all seem to be needed to understand better 
the important problem of crystal growth. 
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